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Abstract 

We propose a theory which describes the density relaxation of loosely packed, 
cohesionless granular material under mechanical tapping. Using the com- 
pactivity concept we develop a formalism of statistical mechanics which allows 
us to calculate the density of a powder as a function of time and compactiv- 
ity. A simple fluctuation-dissipation relation which relates compactivity to 
the amplitude and frequency of a tapping is proposed. Experimental data of 
E.R.Nowak et al. [Powder Technology 94, 79 (1997) ] show how density of 
initially deposited in a fluffy state powder evolves under carefully controlled 
tapping towards a random close packing (RCP) density. Ramping the vibra- 
tion amplitude repeatedly up and back down again reveals the existence of 
reversible and irreversible branches in the response. In the framework of our 
approach the reversible branch (along which the RCP density is obtained) cor- 
responds to the steady state solution of the Fokker-Planck equation whereas 
the irreversible one is represented by a superposition of " excited states" eigen- 
functions. These two regimes of response are analyzed theoretically and a 
qualitative explanation of the hysteresis curve is offered. 



PACS numbers 81.05.R,05.40.,03.20. 
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I. INTRODUCTION 



There is an increasing interest in applying the methods of statistical mechanics and ki- 
netic theory to granular materials where processes are dominated by geometrical constraints 
and friction, and initially posses a memory of sample's preparation f|],@|. 

In this paper we propose an analytical approach which gives a qualitative explanation of 
experimental data obtained by Nowak et al. 0. They have shown that external vibrations 
lead to a slow, approach of the packing density to a final steady-state value. Depending 
on the initial conditions and the magnitude of the vibration acceleration, the system can 
either reversibly move between steady-state densities or can become irreversibly trapped into 
metastable states that is the rate of compaction and the final density depend sensitively on 
the history of vibration intensities that the system experiences (see Fig.l). 

A granular material is a system with a large number of individual grains and therefore it 
has a huge number of degrees of freedom. Grains interact with each other via contact forces 
which are determined by friction, gravitational loading and amplitude of an exernal force if 
the system is perturbed. Therefore one needs to invent a formalism that would allow us to 
calculate macroscopic averages in terms of microscopic (i.e. of individual grains) properties 
of the system. If we assume that it may be characterised by a small number of parameters 
(e.g. analogous to temperature ) and that this system has properties which are reproducible 
given the same set of extensive operations (i.e. operations acting upon the system as a whole 
rather than upon individual grains) then we may apply the ideas of statistical averaging over 
the ensemble of configurations to granular systems M. 

In the present paper we consider the simplest model of a granular material by introducing 
the volume function W and assume the simplest case that all configurations of a given volume 
are equally probable; in many cases the mechanism of deposition will leave a history in the 
configuration but this will not be considered here. W will depend on the coordinates of 
the grains and their orientations and is the analogue of a Hamiltonian. Averaging over all 
the possible configurations of the grains in real space gives us a configurational statistical 
ensemble describing the random packing of grains. Since we are assuming that we are dealing 
with a system whose constituents are hard (i.e. impenetrable) we have to include some 
account of this in our formalism in order to reduce the number of possible configurations 
the system may occupy. Also for a packing which is stable under applied force we must 
consider the configurations restricting the number of possible volume states that the system 
may occupy to be only those configurations which are stable. Also grains cannot overlap 
and this condition produces very strong constraints (frustration) on their relative positions. 
This implies that all grains have to be in contact with their nearest neighbours. Of course 
in the real powder the topological defects can exist such as vacancies, voids or arches. But 
as these will be a subject of a future paper we do not consider them here. Thus we have a 
"microcanonical" probability distribution ||: 

P = e~%5(V - W) Q(contacts) (1) 

e~ = J 5(V — W) ^{contacts) d(all degrees of freedom) (2) 
where we define O as: 
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Q(contacts) 



1 if z > z m 
if z < z m 



where z m is the minimal coordination number of a grain 13. We have to introduce G 
because we consider the stable isotropic and homogeneous packings. Just as in conventional 
statistical mechanics with microcanonical distribution: 



P = e-%6{E-H) (3) 
and temperature: 

T = ^ (4) 

OS y } 

we can define the analogue of temperature as: 

This fundamental parameter is called compactivity HQ. It characterises the packing of a 
granular material and may be interpreted as being characteristic of the number of ways it is 
possible to arrange the grains in the system into volume AV such that the disorder is AS. 
Consequently the two limits of X are and oo, corresponding to the most and least compact 
stable arrangements. This is clearly a valid parameter for sufficiently dense powders because 
one can in principle calculate the configurational entropy of an arrangement of grains and 
therefore derive the compactivity from the basic definition ||. One can expect despite the 
strong constraints resulting from the stability conditions, the number of packings to grow 
exponentially with the volume of a sample and the configurational entropy defined as a 
logarithm of this number is extensive. 

As usual it is more convenient to introduce the canonical probability distribution: 

Y — W 

P = e^r i (6) 

where A is a constant which gives the entropy the dimension of volume, Y we call the effective 
volume, it is the analogue of the free energy: 

e-^=/e-^d(all), V = Y - X^. (7) 

To illustrate this theory consider the simplest example of a W, the analogue of Bragg- 
Williams approximation each grain has neighbours touching it with a certain coordi- 
nation and angular direction. In order to set up an analogy with the statistical mechanics 
of alloys we assume that each grain has a certain property, which defines the "interaction" 
with its nearest neighbours. Taking the coordination number of a grain as such a property 
and assuming that there are just two types of coordination zq and z\ we assign a volume v% 
to any grain with z,i coordination number. Thus we write the volume function as: 

W = n v + (N — noM (8) 
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where N is the number of grains in the system, rii is the number of grains with the 
coordination number Zi and N = no + n\. The simple calculation of Y and V gives us ||: 



Y 



{vq + ^1) 
2 



- iVAX In 2cosh 



AX 



(9) 



V = N 



(vp + vi) 
2 



+ N 



(vq - Vx) 



2 



tanh 



(vq ~ vi) 
XX 



(10) 



Thus we have two limits: V = Nvq, when X — > and V = N(vq + wi)/2 when X —>■ oo 
(X is a number of grains). Note that the maximum V is not Nv\ just as in the thermal 
system (say a spin in a magnetic field) with two energy levels E and E\ one has E = E 
when T — > and E = (Eq + E\)/2 when T — > oo. 

II. "TWO- VOLUMES" MODEL: SOLUTION OF THE FOKKER-PLANCK 

EQUATION 

We consider the rigid grains powder dominated by friction deposited in a container which 
will be shaken or tapped(in order to consider the simplest case we ignore other possible 
interactions e.g. cohesion and do not distinguish between the grain-grain interactions in the 
bulk and those on the boundaries). We assume that most of the particles in the bulk do not 
acquire any non ephemeral kinetic energy i.e. the change of a certain configuration occurs 
due to continuous and cooperative rearrangement of a free volume between the neighbouring 
grains. Any such powder will have a remembered history of deposition and in particular can 
have non-trivial stress patterns, but we will confine the analysis of this paper to systems 
with homogeneous stress which will permit us to ignore it. The fundamental assumption is 
that under shaking a powder can return to a well defined state, independent of its starting 
condition. Thus in the simplest system, a homogeneous powder, the density characterises 
the state. 

It is sensible to seek the simplest algebraic model for our calculation and to this end 
since the orientation of the grain must have at least two degrees of freedom, say \i\ and //2, 
our volume function is: 



implying a two-dimensional picture (see Fig. 2). When fi = we have W = v q then the 
grain is "well oriented" which means that a free volume is minimal and when /i = 1 and 
W = v\ then the grain is "not well oriented" (free volume is maximal). It is a self-consistent 
approximation since the parameters Vq and V\ are the average volumes of the grain in the 
presence of other grains. In general we can write: 



W = v + {v 1 -v )(fj l l + 14) 



(11) 




(12) 



where w(fi) is the weight factor attached to fi. From (|7|) we derive Y and V: 




(13) 
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V = N(v + A.Y) - ^*>. (14) 
e ax — 1 

Thus we have the same limits as for volume function V = iV^o, when X — > and 
V = iV(v + v x )/2 when X -> oo. 

The main physical idea of our approach is the following: all grains in the bulk experience 
the external perturbation as a random force with zero correlation time so that the process 
of compaction can be seen as the Ornstein-Uhlenbeck process for the degrees of freedom 
fj,i, i = 1,2 f7|. Therefore we write the Langevin equation: 

<g + l.™L = jD m (15) 
at v a ^ 

where (fi(t)fj(t')) = 2SijS(t — t') and v characterises the frictional resistance imposed on the 
grain by its nearest neighbours. The term fi(t) on the RHS of fll5|) represents the random 
force generated by a tap. The terms "shaken" or "tapped" have been used above and we 
have to make them more precise. The derivation gives the analogue of the Einstein relation 
that v = (XX)/ D. If we identify / with the amplitude of the force a used in the tapping, 
the natural way to make this dimensionless is to write the "diffusion" coefficient as : 

B=0 2 ^ (16) 
That is we have a simplest guess for a fluctuation-dissipation relation: 

AA ' = (") 2 — (17) 
where v is the volume of a grain, a the frequency of a tap and g the gravitational acceleration. 



Use of the Langevin equation flTT| ) is of course a crude simplification as it does not explicitly 
take into account the presence of boundaries and topological constraints. Generally speaking 
one would have to use the integro-differential Langevin equation with the memory kernel: 



+ f K(t-t') ^(t')dt' = s/Dfi(t), (18) 
Jo 



dt 



as one sees in experiment that the final density depends sensitively on the history of vibration 
intensities. Clearly to solve such an equation is not a trivial task although the solution could 
give us the better understanding of many interesting features of granular compaction. The 
problem of how to choose the initial values of \i is in reality the deposition problem. We 
discuss it later. 

The Langevin equation can be easily solved for W quadratic in /z: 

IH(t) = fJk(0)e-« + VDe-* f f(t')e^' dt' (19) 

J 

Averaging over the ensemble we get: 

(^•(t)> = fH(P)e-* (20) 



5 



where /ij(0) = 1 is the initial value of /ij,7 = 2 ^ Vl ~ v °^ has the meaning of relaxation time 
of the degree of freedom As t —>■ oo fi goes to fif = which corresponds to the random 
close packing limit. The Fokker-Planck equation seems to be quite generic in modelling 
the response of granular materials to an externally applied shear rate || although in that 
problem it is more convenient to use the volume "Hamiltonian" W as a function of the 



coordination number of each grain. The standard treatment of the Langevin equation (15) 
is to use it to derive the Fokker-Planck equation:. 

' Ai^-+T^> = (21) 



dt v dfXidfXj 

where = Ddij and 7^ = 7^. Equation (^1|) can be solved explicitly. It has right- and 
left-hand eigenf unctions P n and Q n and eigenvalues u) n such that: 

d d 

U n Pn = -^—( D iJn 1" lij^j)Pn (22) 

d d 

U n Qn = ( - Dy- h Jijfij) -K—Qn (23) 

or equivalently a Green function: 

G = J EPn(»)QMe- u » t . (24) 

n 

It follows that if we start with a non-equilibrium distribution: 

oo - 

pW(t = 0) = X; A„P„, A n = / Q n p(°) d/i! d/i 2 (25) 

n=0 J 

and it will develop in time as: 

oo 

(t) =A P + J2 A n P n e-" nt (26) 

where J P 1 - ** d/xid//2 = A . This coefficient is determined by a number of grains present in 
the powder, hence must be a constant. The steady-state distribution function is: 

("1~'"0)(^1 +M§) 

p(°) (t oo) = ^ (27) 

Jo e AX d//id^ 2 

The Fokker-Planck operator ([H]) has a complete orthogonal set of eigenfunctions: 

P n = H n e (28) 



where if n are Hermite polynomials and /ij G (0, oo). In our case /ij G (0,1). One 
can avoid this mathematical difficulty, taking into account the crudety of our model and 
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constructing the "first excited state": P 2 = + A 4 ?) + b) e ^ orthogonal to the 
ground-state eigenfunction P . This eigenfunction describes the initial state of our system 
i.e. loosely packed deposited powder. Therefore it is easy to see the initial non-equilibrium 
distribution ( [ZED depends on how the the powder is deposited. Constants a and b can be 
defined from the orthonormality relations. By using: 

P n = QnP m (t^00), Q = l (29) 

and: 

/ Q 2 L FP P 2 dp 1 dp 2 = <jJ 2 (30) 
Jo 

one can easily verify that the eigenvalue u 2 (which corresponds to P 2 and gives us the decay 
rate of our nonequilibrium distribution) is a constant dimensionless number. 

Suppose now that deposition produces a highly improbable configuration, indeed the 
most improbable configuration whith: //? + fi 2 = 2 and the mean volume function is W = 
2t>i — d « Vx, where: 

W(X,t) = J P {0) (X,t)W dprdp 2 . (31) 

It is possible to imagine a state where all the grains are improbably placed, i.e. where 
each grain has its maximum volume v\. In a thermal analogy this would be like fully 
magnetised magnetic array of spins where the magnetic field is suddenly reversed. Such a 
system is highly unstable and equilibrium statistical mechanics does not cover this case at 
all. It will thermalize consuming the very high energy whilst establishing the appropriate 
temperature. Powders however are dominated by friction, so if one could put together a 
powder where the grains were placed in high volume configuration, it will just sit there until 
shaken; when shaken it will find its way to the distribution (|BJ). It is possible to identify 
physical states of the powder with characteristic values of volume in our model. The value 
V = Nvi corresponds to the "deposited" powder, i.e. the powder is put into the most 
unstable condition possible, but friction holds it. When V = Nvq the powder is shaken into 
closest packing. The intermediate value of V = (v + i>i)/2 corresponds to the minimum 
density of the reversible curve. Thus we can offer an interpretation of three values of density 
presented in the experimental data 0. 

The general solution of the Fokker-Planck equation fl2~l|) goes to its steady-state value 
when t — > oo so we can expect W(X,t) to diminish (as the amplitude of tapping increases) 
until one reaches the steady-state value W(X). The formula ([14] ) can be obtained using fl3~l]) 
when t — > oo and represents a reversible curve in experimental data of ||: altering a moves 
one along the curve p = y/(x) = P( a )- We can identify time with the number of taps, so 
wherever we start with any initial p(o) and a, successive tapping takes one to reversible curve 
p(a). Or, if one decides on a certain number of taps, t ^ oo , one will traverse a curve pt(a), 
where Poo(a) = p{o). Notice that the simple result lies within the crudety of our model. The 
general problem will not allow us to think of X as X(a) independent of the development of 
the system. The thermal analogy is this: if the Brownian motion in an ensemble of particles 
is controlled by a random force / which is defined in terms of its amplitude and time profile, 
this random force defines the temperature in the system. Our problem is like a magnetic 
system where magnetic dipoles are affected by a constant magnetic field, being random at 
high temperature , and increasingly oriented by the external field as the temperature falls. 
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III. DISCUSSION 



The physical picture presented in section 2 is consistent with everyday knowledge of 
granular materials: when poured they take up a low density but when shaken settle down, 
unless shaken violently when they return to low density. These effects are much more pro- 
nounced in systems with irregularly shaped grains then with fairly smooth uniform spheres, 
indeed the more irregular a grain is, the more the discussion above describes big differ- 
ences between p( ) and p. The experimental data of || show the packing density depen- 
dence on parameter V = a/g for a fixed number of taps. A loosely packed bead assembly 
first undergoes irreversible compaction corresponding to the lower branch of p{T). The 
settling behavior becomes reversible only once a characteristic acceleration has been ex- 
ceeded. Our theory gives three points p(X = 0),p(X = oo) and p{t = 0) which are in 
the ratio: Vq 1 , , ^ > , ff 1 and these are in reasonable agreement with experimental data: 

p(X = 0) = pa 0.64, p = ^ ~ 0.58 and p(X = oo) = (^p^y ~ 0.62. Another important 
issue is the validity of the compactivity concept for a "fluffy" but still mechanically stable 
granular arrays e.g. for those composed of spheres with p < 0.58. In our theory p(X = oo) 
corresponds to the beginning of the reversible branch (see Fig.l) and using our analogy with 
a magnetic system is analogous to dipoles at a high temperature. The irreversible branch 
has an analogue in the behaviour of the magnetic system where initially the dipoles are 
strongly aligned with an external field but this field is then flipped to the opposite direction. 

The fluffy powder is a very complicated object as it has plenty of topological defects 
and stress arches. Throughout the paper we assumed that our granular array is spatially 
homogeneous which is the case for densities of the reversible curve. However this is a very 
subtle problem which will be a subject of a future paper. It is a difficult problem to decide 
whether embarking on a vast amount of algebraic work that a superior mode would entail is 
worthwhile. But our simple model is quite physical and can be extended when experiments 
would justify the work. 

A final point is that we find the lower (irreversible) curve build up to the upper (re- 
versible) curve exponentially in time: 

V{t) = initial e-* + Wl - e" *) (32) 

while one can expect the logarithmic in time approach to the steady state density e.g. 
the Vogel-Fulcher type curve which is typical of disordered thermal systems such as spin 
and structural glasses || 

V(t) = V f + (Vi - Vf)e- Ut + Tt~ € (33) 

where e is large. 

In fact, identifying t with the number of taps n, the law seems to be even slower at (lnt) -1 . 
Our simple analysis is clearly inadequate to obtain such a result which is quite outside the 
straightforward method of expansion in the present set of eigenfunctions. However there is 



an argument by de Gennes |10| which argues that a Poisson distribution can provide this 
logarithmic behaviour. 
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FIGURES 




FIG. 1. Dependence of the steady-state packing density on the tapping history (Nowak et al.). 
Experimental values of density packing fraction are in the following correspondence with model 
parameters: p(X = 0) = ^ « 0.64, p(t = 0) = p = ± « 0.58 and p(X = oo) = w 0.62. 





(1=0 



FIG. 2. Graphical representation for the limit values of the degree of freedom p in 2-D 
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